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一 种 改进 的 警示 传播 算法 求解 Max-SAT 问题 
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摘 要 : Max-SAT 问题 是 SAT 问题 的 优化 版 本 ,目标 是 在 给 定 的 子 句 集中 ,找到 一 组 变 元 赋值 ， 使 得 满足 子 句 数 最 
多 ， 该 问题 是 典型 的 NP-hard 问题 。 随 着 大 数据 和 人 工 智能 的 深度 发 展 ， 过 去 原 有 的 算法 已 不 再 适用 ,设计 新 的 求 
解 算法 或 对 已 有 的 求解 算法 进行 优化 是 目前 研究 的 热点 。 本 文 针 对 警示 传播 算法 求解 随机 Max-3-SAT 问题 的 局 限 
性 ， 提 出 了 一 种 基于 变 元 权 值 计 算 的 警示 传播 算法， 结合 随机 游 走 算法 ， 给 出 一 种 新 型 算法 WWP+WalkSAT， 通 过 
改进 求解 的 局 限 性 ， 更 好 地 得 到 一 组 有 效 的 初始 解 ， 从 而 提高 算法 的 局 部 搜索 能 力 。 利 用 2016 年 Max-SAT 国际 竞 
赛 部 分 基准 实例 ， 将 WWP+WalkSAT 算法 与 8 种 局 部 搜索 算法 进行 精度 方面 的 对 比 实验 。 实 验 结 果 表 明 
WWP+WalkSAT 算法 有 较 好 的 性 能 。 

关键 词 : 可 满足 性 问题 ; 最 大 可 满足 性 问题 ; 警示 传播 算法 ; 局 部 搜索 算法 

中 图 分 类 号 : TP301 doi: 10.19734/j.issn.1001-3695.2022.01.0023 


Improved warning propagation algorithm for solving max-sat problem 


Wu Yuxiang’, Wang Xiaofeng*", Ding Hongshengsi , Yu Zhuos 
(a. School of Computer Science & Engineering, b. The Key Laboratory of Images & Graphics Intelligent Processing of State 
Ethnic Affairs Commission, North Minzu University, Yinchuan 750021, China) 


Abstract: The Max-SAT problem is an optimized version of the SAT problem. The goal is to find a set of variable assignments 


in a given set of clauses so that the maximum number of clauses is satisfied. This problem is a typical NP-hard problem. With 
the in-depth development of big data and artificial intelligence, the original algorithms in the past are no longer applicable. 
Designing a new solution algorithm or optimizing the existing solution algorithm is the current research hotspot. Aiming at 
the limitations of the information dissemination algorithm for solving the random Max-3-SAT problem, this paper proposes a 


warning dissemination algorithm based on the calculation of variable weights. Combined with the random walk algorithm, a 
a bs new algorithm WWP+WalkSAT is formed, which is solved by improvement The limitation of, it is better to get a set of 
effective initial solutions, thereby improving the local search ability of the algorithm. Using some benchmark examples of the 
re Max-SAT International Competition in 2016, the WWP+WalkSAT algorithm and 8 local search algorithms are used for 
accuracy comparison experiments. The experimental results show that the WWP + walksat algorithm has good performance. 


Key words: satisfiability problem; max-sat problem; warning propagation algorithm; local search algorithm 


全 0 引言 主要 分 为 完备 性 算法 和 非 完备 性 算法 两 种 。 完 备 性 算法 是 可 
二 以 找到 问题 的 精确 解 ， 但 在 求解 大 规模 问题 时 ， 由 于 其 复杂 

组 合 优 化 问题 在 运筹 学 、 离 散 数 学 和 计算 机 科学 中 都 有 度 是 指数 级 ， 很 难 在 合理 时 间 找 到 结果 。 近 年 来 ， 人 们 主要 
着 非常 重要 的 作用 ， 并 在 国防 、 交 通 、 医 疗 、 通 信 等 领域 有 在 分 支 策略 、 推 理 规则 、 下 界 估计 世 5 三 方面 进行 改进 ， 进 
着 广泛 的 应 用 ， 其 中 常见 的 组 合 优化 问题 包括 背包 问题 而 出 现 了 许多 有 效 的 算法 ， 具 有 代表 性 的 有 : WmaxSatz[0、 
(Knapsack)、 旅 行商 问题 (Travelling salesman problem, TSP)、 MiniMaxSat03 等 。 相 对 于 完备 算法 而 言 ， 非 完备 性 算法 则 是 
车 辆 路 径 规划 问题 (Vehicle Routing Problem,VRP)、 最 大 团 问 能 短 时 间 内 找到 大 规模 问题 的 最 优 解 , 从 而 提高 了 求解 效率 ， 
题 (Maximum Clique Problem，MCP) 、 最 小 顶点 覆盖 问题 但 缺点 则 是 不 能 保证 解 的 准确 性 。 主 流 的 非 完 备 性 算法 是 近 
(Minimum vertex cover MVC)、 最 大 可 满足 性 问题 (Maximum 以 算法 、 信 息 传播 算法 、 局 部 搜索 算法 03-5]。 
Satisfiability Max-SAT0-9 等 。 而 最 大 可 满足 性 问题 是 一 个 信息 传播 算法 是 一 种 启发 式 信息 传递 算法 ， 起 源 于 统计 


型 的 NP 难 问 题 , 是 可 满足 性 问题 (Satisfiability， SAT)JI 的 优 ”物理 学 ， 该 类 算法 通过 边缘 概率 计算 的 方法 ， 已 应 用 在 许多 
化 版 本 。 给 定 一 个 命题 合 取 范式 (CNF)， CNF 是 由 一 组 子 句 “领域 。 特 别 是 对 于 命题 公式 的 可 满足 性 问题 ， 有 三 种 基于 因 
合 取 构 成 ， 每 一 个 子 句 由 一 组 变 元 析 取 构成 ，Max-SAT 问题 。” 子 图 的 信息 传播 算法 ， 即 警示 传播 算法 (warning propagation， 
是 寻找 一 组 赋值 使 得 满足 子 句 数目 最 多 。Max-SAT 在 组 合拍 WP)、 信 念 传播 算法 (belief propagation, BP) 和 调查 传播 算法 
卖 、 车 辆 调度 、 排 课 安排 等 现实 问题 中 有 着 广泛 的 应 用 ， 后 (survey propagation, SP)。 目 前 ， 在 求解 随机 SAT 问题 上 信息 
时 图 论 中 的 最 大 团 问 题 和 顶点 支配 集 也 可 以 转换 成 Max- 传播 算法 是 最 为 肥效 的 方法 , 并 且 在 图 着 色 、 最 大 流 、LDPC 
SAT 问题 进行 求解 。 编码 等 许多 组 合 优化 问题 中 取得 了 不 错 的 成 果 。 

为 了 解决 Max-SAT 问题 ， 研 究 人 员 提 出 许多 有 效 算法 ， 然而 ， 信 息 传播 算法 的 收敛 性 和 有 效 性 一 直 是 研究 的 重 
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点 。 收 敛 性 是 


言 恩 传播 算法 两 次 的 信息 迭代 值 不 发 生变 
化 ， 有 效 性 是 指 : 信息 传播 算法 能 够 有 效 的 求解 问题 。 文 献 
[16~18] 分 析 了 信息 传播 算法 的 收敛 性 和 有 效 性 ,并 给 出 了 算法 
改 敛 的 充分 条 件 。 文 献 [19] 则 是 基于 具体 的 实例 产生 模型 分 析 了 
言 恩 传播 算法 的 收敛 性 ,并 给 出 了 算法 收敛 的 概率 条 件 。 
进一步 研究 发 现 , 在 随机 3-SAT 问题 中 , 子 句 个 数 m 和 
变 元 个 数 n 的 比值 称 约束 密度 & , 即 a=m/n，, 它 不 仅 对 公式 
的 可 满足 性 产生 影响 ， 还 对 可 满足 性 公式 的 求解 难 易 程度 有 
着 非常 重要 的 影响 。 随 着 约束 密度 的 增 大 ,， 当 3.52<w<448 
时 ， 会 发 生 相 变 。 在 相 变 范围 以 外 的 可 满足 性 实例 均 为 易 解 
实例 ， 高 概率 是 可 满足 的 ;在 相 变 点 附近 的 的 实例 属于 难 解 
实例 ， 高 概率 是 不 可 满足 的 。 虽 然 信 息 传 播 算法 对 于 求解 难 
解 实例 十 分 有 效 ， 但 在 求解 相 变 点 以 外 的 易 解 实例 方面 ， 信 
息 传播 算法 不 能 收敛 。 

针对 于 此 问题 ， 文 本 设计 了 一 种 新 型 变 元 权 值 计 算 的 警 
示 传 播 算法 WWP+WalkSAT， 通 过 变 元 权 值 计算 对 警示 传播 
算法 进行 改进 ， 可 以 得 到 一 组 有 效 的 初始 解 ， 再 结合 随机 游 
走 算法 进行 局 部 搜索 。 对 该 算法 的 性 能 表现 进行 实验 分 析 ， 
在 相 变 点 附近 实例 和 相 变 点 以 外 实例 均 与 其 余 局 部 搜索 算法 
进行 实验 对 比 。 实 验 结果 表明 : 对 于 求解 任何 区 域 的 实例 ， 
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WWP+WalkSAT 在 求解 精度 方面 平均 优 于 其 他 局 部 搜索 算法 。 


1 ”基础 知识 


最 大 可 满足 性 问题 是 指 : 给 定 一 组 命题 变 元 ssX ,由 这 
些 变 元 形成 一 组 子 句 ， 构 成 CNF 公式 ， 使 得 CNF 公式 中 满 


足 子 句 的 个 数 最 多 ， 即 ， 使 得 不 满足 的 子 句 个 数 最 少 。Max- 
SAT 问题 的 数学 模型 如 式 (1) 和 式 (2): 
Maximize: 
2 “Za (1) 
Subject to: vaeF: +t) 
z, Ee{0,1} x e{0,1} (2) 


文字 的 集合 (Literal) 关 ={%,w, 累 ,…, 加 }: 每 一 个 布尔 变 
元 neX ， 变 元 # 取 正 时 表示 正文 字 ; 变 元 * 取 反 时 表示 负 
文字 。 

子 句 的 集合 (Clause) C={c,C,C…,Cj : 对 m 个 不 同 的 子 
名 形成 一 个 子 名 集合， 每 一 个 子 句 由 一 个 或 多 个 文字 组 成 ， 
子 句 与 子 句 之 间 通 过 析 取 运算 连接 ， 且 仅 当 子 句 中 至 少 有 一 
个 文字 为 1 时， 该 子 句 满足 ， 反 之 子 句 不 满足 。 每 一 个 子 句 
中 文字 的 个 数 叫做 这 个 子 名 的 长 度 。 当 子 句 中 只 有 一 个 文字 
出 现时 ， 称 该 子 名 为 单子 句 。 

合 取 范式 (CNF 公式 ): 由 若干 个 子 句 合 取 构 成 ， 即 
G 人 CG 和 人 G 人-… 人 和 ^G，， 当 且 仅 当 CNF 公式 的 每 一 个 子 句 都 满足 
时 ， 称 合 取 范式 CNF 可 满足 。 

因子 图 (Factor Graph): 因子 图 也 称 二 部 图 ， 图 中 包含 两 
类 节点 : 一 类 是 变 元 节点 (图 中 用 圆 环 表 示 ， 记 住 (1，2，3， 
4;…)), 男 一 类 是 函数 节点 (图 中 中 方 框 表示 , 记 住 (a,b,c;*…))。 
给 一 个 CNF 公式 f={G,G,G,.…G,} ,公式 中 由 ,加 ,加 ,…h 个 变 
元 构成 ， 其 中 所 有 变 元 对 应 的 是 变 元 节点 ， 变 元 析 取 构成 的 
子 句 对 应 的 是 函数 节点 。 图 中 出 现 的 虚 边 和 实 边 代 表 着 CNF 
公式 的 两 种 不 同 的 含义 : 实 边 : 子 句 中 代表 正文 字 罗 , 虚 边 : 
子 句 中 代表 负 文 字 5 。 

例 1: F={Wvwv A VE VA MVGVG)ARG Vv))}= 
taApAcAd} 其 对 应 的 因子 图 如 图 1 所 示 。 

VD=V*(W)+V-(@) 表示 子 句 a 中 出 现 的 变 元 集合 
V(@)=Vi(@)+V-(a) 。 其 中 : V*(@) 表示 子 句 a 中 变 元 为 正 的 集 
合 ，Y-(O 表示 子 句 a 中 变 元 为 负 的 集合 ，VCOV 表 示 子 句 a 
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中 除 变 元 i 以 外 其 余 变 元 的 集合 ; 在 图 1 中 ，Y(O ={amz} 
Vi(W={0)}, VW)={%,%} 

V() 表示 含有 变 元 怀 的 子 句 集合 YCOD=CD+Y CD， 其 中 
CD 表示 变 元 沁 为 正 的 子 句 集合 ，Y-(0D 表示 变 元 5 为 负 的 


子 句 集合 ，VYCh\a 表示 除了 子 句 a 以 为 ， 其 余 含 有 变 元 5 的 
子 句 集合 。 在 图 1 中 ，YO ={wpecl ,V+(D={a},V-(D)={b,c} 。 


图 1 因子 图 


Fig. 1 Factor graph 
及 :对 于 子 句 a 中 的 变 元 5 的 标识 符 , 变 元 在 子 句 中 的 取 值 ， 
由 式 (3) 表 示 为 


一 1 XEa 
| —xjea G3) 

定义 两 个 一 致 性 集合 WG) 和 Vi()) 

W() 表示 除 子 句 a 以 外 , 含有 变 元 且 变 元 二 的 取 值 方 
向 与 其 子 句 a 中 取 值 不 一 致 的 子 句 集 ; 让 (CD 表示 除 子 句 a 以 
外 ， 含 有 变 元 罗 且 变 元 ;的 取 值 方向 与 其 子 句 a 中 取 值 一 致 
的 子 句 集 ; 集合 WW()) 和 到 (六 如 图 2 所 示 。 


Vi (k) Ve (k) 
/m9 I 
学 
Vy 
k ¥-- 


图 2 局 部 因子 图 
Fig.2 Part of factor graph 


2 ”基于 变 元 权 值 计算 的 警示 传播 算法 


基于 变 元 权 值 计算 的 警示 传播 算法 WWP+WalkSAT 是 
在 警示 传播 算法 PCO 上 进行 改进 ， 文 献 [16] 分 析 警 示 传 播 算法 
的 收敛 性 ， 通 过 引入 参数 ， 将 算法 中 的 信息 取 值 从 {0.1 松弛 
为 [0 ， 利 用 向 量 空间 中 压缩 函数 的 性 质 ， 给 出 警示 传播 算 
法 收敛 的 一 个 充分 条 件 ， 并 为 警示 传播 算法 性 能 提供 了 理论 
依据 。WWP+WalkSAT 是 在 此 理论 基础 上 ， 引 入 了 变 元 权 值 
计算 Ri， 从 中 得 到 一 组 有 效 的 初始 解 ， 其 目的 是 选 出 权 值 最 
高 的 文字 ， 使 得 满足 子 句 个 数 最 多 ， 从 而 减少 无 效 的 子 句 传 
播 ， 加 快 搜索 过 程 。 原 本 的 警示 传播 算法 在 求解 难 解 实例 十 
分 有 效 ， 但 对 于 易 解 实例 会 出 现 不 收敛 现象 。 对 于 此 问题 ， 
该 算法 在 迭代 次 数 结束 时 ， 计 算 局 部 警示 信息 变化 最 小 时 的 
警示 值 ， 再 进行 变 元 权 值 计算 。 最 后 对 其 初始 解 采用 随机 游 
走 的 局 部 搜索 算法 ， 可 得 到 更 优 解 ， 提 从 而 高 其 算法 性 能 。 
WWP+WalkSAT 打破 了 原先 警示 传播 算法 求解 Max-SAT 问 
题 的 局 限 性 ， 可 用 于 求解 任何 区 域 的 实例 ， 同 时 在 求解 精度 
上 也 有 所 提高 。 
2.1 警示 传播 算法 (Warning Propagation) 

通过 信息 传递 而 设计 的 信息 传播 算法 ， 对 于 求解 可 满足 
性 问题 具有 良好 的 有 效 性 ， 警 示 传 播 算法 是 一 个 迭代 运算 的 
算法 ， 在 每 一 次 迭代 过 程 中 ， 对 于 因子 图 的 每 条 边 (a,? ， 都 
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会 获得 一 个 警示 信息 ， 表 示 在 因子 图 中 一 个 函数 节点 a 传递 
给 变 元 i 一 个 布尔 信息 , 记 作 ww sf01 。w> 的 迭代 方程 如 式 
(4) 所 示 。 
C0) (4) 
JEV(aJV beV(jN\a 


时 中 f 是 迭代 次 数 ，eco 是 一 个 截 尾 函 数 ， 定 义 为 : 
oa 二 9， 由 此 公式 可 以 说 明 ， 当 w=-1 表 示 对 于 子 句 


—l], x<0 
a 的 满足 性 由 变 元 所 来 决定 ， 当 ww =0 表 示 子 句 a 的 满足 性 
不 能 完全 由 变 元 二 决 定 ， 即 子 句 a 的 满足 性 由 其 余 变 元 的 取 
值 来 决定 。 
对 于 每 一 个 变 元 二， 需要 计算 一 个 局 部 腔 域 也 和 一 个 冲 


突 域 4e{0,1}， 其 中 局 部 腔 域 和 冲突 域 式 (5) 和 式 (6) 所 示 。 
ee J ? byi = pHi— bi 


lif( 2 wi)( 2 wi)>0. 
Ci= beVi*(i) beV-(i) (6) 


0, otherwise. 


当 冲 突 域 4=1 时 ， 即 4, 人 多 ,“)>0 ， 则 表示 对 于 


变 元 ,受到 了 两 个 子 句 的 约束 ， 变 元 在 两 个 子 句 中 的 取 
值 不 一 致 ， 因 此 变 元 * 的 取 值 发 生 了 冲突 。 如 : 子 句 a 告诉 
变 元 所 取 1， 子 句 b5 告诉 变 元 入 取 0。 当 冲突 域 oc=0 时 ， 即 


(之 发 2 区 )=0，, 则 表示 对 于 变 元 的 取 值 没 有 冲突 发 生 ， 


beV+ (i) beV-(i) 


变 元 的 具体 取 值 可 以 通过 局 部 腔 域 来 决定 ， 如 式 (7) 所 示 。 
1 > 


beV+(i) beV-(i) 


5=1(00 > wy = Dw,=0 (7) 
beV+ (i) beV- (i) 
0, > u;,; =0, 2 wy >0 


2 bev) 
如 果 五 >0， 变 元 =1; 如 果 矿 <0， 变 元 *=0; 否则 变 
元 * 暂时 不 赋值 。 
算法 1 警示 传播 算法 (WP) 
输入 : CNF 因子 图 6; 最 大 迭代 次 数 tx ; 精度 2 。 


输出 : 未 收敛 或 收敛 ， 输 出 警示 信息 值 zs 
Begin 
当 !=0 时 ， 对 于 因子 图 的 每 个 边 (a,i) ， 以 1/2 的 概率 随机 初始 化 
警示 信息 w_,(0) 一 {01 // 随 机 初始 化 
For t =1 to ta 
通过 字典 的 顺序 ， 对 图 中 每 一 条 边 (a,i) eE 进行 遍历 ， 采 用 警示 

更 新 算法 UPDATE-WP 对 w， 即 w(t—D uD) 


; 进行 更 新 ， 


www(-DI<e ， 说 明 


当 图 中 所 有 的 边 (4,i)eE， 若 


已 收 仇 ， Wy(D=wWy(t-D， 
LIF f= 


停止 循环 ， 返 回 ww(7) 


寻找 六 =argmin 2 |usw(D) -uo y(t-D)| 


1 (aieE 


返回 : 未 收敛 ,并 将 
并 赋值 给 警示 值 
算法 2 ”警示 信 
输入 : 警示 信息 值 集合 {wj:beV(D\a }， 
输出 : 新 的 wy 


// 把 x* 时 刻 的 信息 值 保 


ai(D = u(t) , 


现 


息 值 更 新 算法 (UPDATE-WP) 
对 所 有 的 jeV(@D)\i 


Begin 
IF: V(D\i=J ， 则 返回 ,=1 
ELSE: 
for jeV(D\i 
IF V()\a =G, 则 hs,=0 
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Else hy = (Boor, cpuls ss; )— (Tour rat si) 
返回 wy =Tjoyew0( Yh ) 
算法 3 曾 示 信息 赋值 指派 和 法 
输入 : CNF 因子 图 G 
输出 : 输出 一 组 变 元 赋值 指派 {7, 怠 , 吉 … 克 ,} 
Begin 
运行 WP 算法 
计算 所 有 的 局 部 腔 域 8; 和 冲突 域 c 
IF c=0， 则 对 变 元 进行 权 值 计 算 


IF c=1 
IF H;>0, 则 x =1 
IF HH,<0, 
ELSE Hi;=0 , 则 对 变 元 进行 权 值 计算 
返回 变 元 赋值 指派 
2.2 变 元 权 值 计算 
在 迭代 运算 时 ， 基 础 的 警示 传播 算法 只 对 冲突 域 " = 

且 局 部 腔 域 五 >0 的 变 元 赋值 为 1， 局 部 腔 域 厂 <0 ee 
直 为 0, 而 对 于 冲突 域 6=1 的 变 元 和 c=0 且 =0 的 变 元 无 法 
起 值 。 在 求解 最 大 可 满足 性 问题 中 ， 变 元 赋值 的 好 坏 直 接 影 
响 着 整个 算法 的 求解 效率 ， 正 确 的 变 元 赋值 可 以 高 效 求解 。 
本 文 提出 了 一 种 新 的 启发 式 变 元 权 值 计 算 ， 通 过 权 值 计算 来 
确定 该 变 元 的 赋值 。 
对 于 未 赋值 的 每 一 个 变 元 都 有 一 个 正 或 负 的 权 值 ， 这 取 
决 于 该 变 元 的 正文 字 和 人 负 文 字 出 现 的 差异 ， 变 元 的 权 值 代表 
着 变 元 为 正 或 为 负 的 程度 。 如 果 权 值 为 正 ， 则 代表 变 元 在 整 
个 CNF 公式 中 更 多 地 以 正文 字 的 形式 出 现 , 而 不 是 负 文 字 的 
形式 出 现 ， 反 之 亦 然 。 权 重 计算 如 式 (8): 
Number( PosLit) — Number( NegLit) (8) 

Number(Clause,,, ) 


其 中 Numper(PosLit) 表示 变 元 正文 字 出 现 次 数 ， Number(NegLit) 
表示 变 元 负 文 字 出 现 次 数 ， Number(Clausew) 表示 变 元 在 子 句 
中 出 现 的 次 数 。 

算法 4 ” 变 元 权 值 算法 
输入 : CNF 公式 

输出 : 变 元 x 的 赋值 


则 =0 


I 


受 寥 


Wer 


[ 


For 对 每 一 个 未 赋值 的 变 元 i: 
初始 变 元 权 值 W =0 ， 找 到 包含 变 元 i 的 所 有 子 句 Clause; 并 记 
录 总 子 句 数 Number; 
For 对 每 一 个 子 句 Clause; 中 的 变 元 
If i>0, 
Elsei<0, 


PosLit ++ 
NegLit ++ 
PosLit — NegLit 
Number 
则 返回 闷 =1 
Else 返回 x=0 
2.3 局 部 搜索 算法 
在 局 部 搜索 算法 中 ， 搜 索 过 程 占用 整个 算法 的 大 部 分 
时 间 ， 普 遍 随机 搜索 算法 的 初始 解 都 是 随机 选取 ， 会 造成 
初始 解 的 随机 性 ， 好 的 初始 解 可 以 高 效 的 找到 结果 ， 不 好 
的 初始 解 会 浪费 很 长 时 间 且 降低 了 算法 的 性 能 ， 同 时 还 会 
重复 地 访问 一 些 解 ， 出 现 循环 的 现象 。 本 文 提出 的 变 元 权 
值 的 警示 传播 算法 ， 可 以 有 效 得 到 一 组 初始 解 ， 对 后 续 进 
行 局 部 搜索 有 很 大 帮助 。 
于 变 元 权 值 的 警示 传播 算法 得 到 一 组 变 元 的 赋值 ， 
此 构造 了 一 组 初始 解 。 初 始 解 得 到 的 子 句 数 不 一 定 是 最 优 的 ， 
需要 局 部 搜索 算法 进行 深 一 步 搜索 。 对 于 Max-SAT 问题 ， 
的 是 找到 满足 子 句 数 最 多 的 一 组 赋值 ， 将 得 到 的 一 组 变 元 赋 


W= 


IF W>0, 


立 | 
口 


虹 
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值 放 入 到 WalkSAT 中 进行 搜索 。 法 对 于 求解 易 解 算 例 十 分 有 效 ， 而 且 优 于 普通 的 局 部 搜索 算 
法 。 表 3 中 ,实验 采用 HG-3SAT-V300-C1000、HG-3SAT-V250- 


Sh 
3 ”实验 结果 与 分 析 C1200 两 种 数据 集 (3SAT 实例 ，300 个 变 元 ，1000 和 1200 个 
为 了 准确 评估 本 文 提出 的 新 型 变 元 权 值 计算 的 警示 传播 算 。 子 句 ，a=4) 。 实验 结果 表明 , 算法 CCLS2akms 是 无 法 求解 
法 能 在 解决 Max-SAT 问题 时 起 到 很 好 的 效果 ， 本 文 对 算法 该 类 问题 的 ，WWP+WalkSAT 的 求解 精度 均 达 到 了 98~99% 
WWP+WalkSAT、 WalkSAT、 SA、 GA、 VNS-GA、 CCLS2akms、 左右 ， 而 CCEHC 的 精度 稳定 在 99% 左 右 。WWP+WalkSAT 
CCEHC、Optirise6-in 和 HS-Greedy 进行 了 测试 ,其 中 WalkSAT、 还 是 略 差异 CCEHC， 但 是 性 能 远 远 好 于 SA、WalkSAT、 
SA、GA、VNS-GA 为 经 典 的 启发 式 算法 ,CCLS2-akms、CCEHC、 ”Optirise6-in 和 HS-Greedy 几 种 算法 。 在 难 解 算 例 中 ， 
Optirise6-in 和 HS-Greedy 为 2016 年 Max-SAT Evalution 国际 竞 WWP+WalkSAT 稍 差 于 CCEHC 算法 ， 但 在 易 解 实例 中 ， 
赛 的 求解 器 。 实 验 结果 使 用 了 2016 年 Max-SAT Evalution 国际 WWP+WalkSAT 求解 性 能 好 于 CCEHC 算法 。 由 此 可 以 说 明 
竞赛 中 随机 类 的 数据 集 。 WWP+WalkSAT 算法 对 求解 随机 Max-3-SAT 的 问题 十 分 有 交 
表 1 给 出 了 SA、WalkSAT、WWP+WalkSAT、GA、VNS- 除 此 之 外 ， 通 过 5 个 数据 集 对 算法 WWP+WalkSAT 和 有 
GA 几 种 算法 之 间 的 统计 比较 结果 。 采 用 s3v70c1000、 WalkSAT 在 迭代 次 数 方面 进行 了 对 比 实验 。 图 3~7 表明 ,7 
s3v80c1000 两 种 数据 集 (3-SAT 实例 ， 变 元 个 数 70 和 80,， 子 者 均 进行 1000 次 迭代 , 结果 表明 除了 数据 集 S3v110c1000 中 
名 个 数 1000), 对 满足 最 大 子 句 个 数 的 结果 进行 了 对 比 实验 。 初始 结果 WWP+WalkSAT 略 差 于 WalkSAT, 其 余 四 个 数据 集 
表 中 可 以 看 出 ， 在 70 个 子 句 中 的 10 个 实例 中 ， 有 7 个 实例 ”一 开始 WWP+WalkSAT 就 优 于 WalkSAT。 这 是 因为 WalkSAT 
都 是 取得 最 好 的 结果 ， 有 3 个 实例 略 差 一 些 。 在 80 个 子 名 的 初始 解 是 随机 赋值 ， 当 初始 解 赋值 好 的 时 候 会 优 于 
中 ， 只 有 4 个 实例 取得 最 好 结果 ， 其 余 6 个 效果 稍 差 一 些 。 WWP+WalkSAT, 而 WWP+WalkSAT 的 初始 解 是 通过 变 元 权 
但 是 可 以 明显 看 出 WWP+WalkSAT 的 效果 远 远 高 于 1 


lu 


a 


直 计 算得 出 的 ， 所 以 一 开始 的 满足 子 句 个 数 就 会 优 于 
WalkSAT， 且 满足 的 子 句 个 数 也 大 大 提高 。 WalkSAT。 而 随 着 迭代 次 数 的 增加 ，WWP+WalkSAT 满足 子 
表 1 SA、WalkSAT、WWP+WalkSAT、GA 和 名 的 个 数 远 远 优 于 WalkSAT。 尽 管 在 数据 集 S3v90c900 中 出 
VNS-GA 在 易 解 实例 中 实例 对 比 现 WalkSAT 高 于 WWP+WalkSAT 的 现象 ， 但 是 最 后 
Tab.1 ”Comparison of SA, WalkSAT, WWP+WalkSAT,GA and WWP+WalkSAT 还 是 取得 了 不 错 的 结果 。 虽 然 在 精度 上 取得 
VNS-GA in easily solved instances 了 不 错 的 效果 , 但 是 WWP+WalkSAT 算法 的 运行 时 间 会 差 于 
The number of satisfy clauses WalkSAT 算法 ， 这 也 是 该 算法 今后 需要 改进 的 地 方 。 
nn SA WalkSAT WWP+WalkSAT GA VNS-GA 福 td a 大 hd - 
s3v70c1000-1 ”950 929 950 921 946 § | $ A 
s3v70c1000-2 948 924 955 931 952 : 
s3v70c1000-3 945 927 950 922 949 oo lb 
s3v70c1000-4 948 930 948 928 948 和 生 
s3v70c1000-5 951 926 951 933 955 = | E 
s3v70c1000-6 ”946 926 948 929 945 EE 
s3v70c1000-7 955 929 953 930 943 | Ta ee 
s3v70c1000-8 947 931 953 931 950 图 3 数据 集 : s3v70c700 图 4 数据 集 : s3v90c900 
s3v70c1000-9 955 925 956 930 958 Fig.3 Dataset: s3v70c700 Fig.4 Dataset: s3v90c900 
s3v70c1000-10 949 930 954 928 952 tt 人 i ee 
s3v80c1000-1 。 953 926 955 953 953 $ er) = 
s3v80c1000-2 951 925 942 953 953 2 ; 
s3v80c1000-3 958 930 949 946 955 : < 
s3v80c1000-4 942 929 953 951 950 3 
s3v80c1000-5 954 927 953 955 955 : 2 
s3v80c1000-6 950 934 958 956 956 本 本 i 
s3v80c1000-7 948 927 953 957 957 ea aa 
s3v80c1000-8 952 929 950 955 955 图 5 数据 集 : s3v110c1000 ”图 6 数据 集 : HG-3SAT-V250-C1000 
s3v80c1000-9 950 932 953 958 959 Fig. 5 Dataset: S3v110c1000 Fig.6 Dataset HG-3SAT-V250-C1000 
s3v80c1000-10 ”947 931 959 952 957 mm DO 
表 2~3 是 将 SA、WalkSAT、WWP+WalkSAT、CCEHC、 : 
CCLS2akms、Optirise6-in 和 HS-Greedy 进行 了 测试 。 实验 结 a | 
果 使 用 了 2016 年 Max-SAT Evalution 竞赛 中 随机 类 的 数据 集 。 | 
其 中 表 2 中 采用 3-SAT 实例 ， 变 元 数 70、90、110 个 ， 子 名 六 
数 700、800、900、1000、1100 的 数据 集 。 每 一 种 实例 有 10 a™"| 3 
组 数据 。 实 验 结果 采用 平均 精度 计算 表示 ， 即 最 大 可 满足 子 i I 
句 数 / 实 例子 句 数 ， 再 对 10 组 数据 求 平 均值 。 实 验 表 明 : 对 0 
于 WWP+WalkSAT 来 说 , 精度 达到 了 95%~99%, 而 CCEHC 图 7 数据 集 : HG-3SAT-V300-C1200 
的 精度 在 95% 到 97%。 且 在 15 种 数据 集中 ， 有 12 个 数据 集 Fig.7 Dataset: HG-3SAT-V300-C1200 


WWP+WalkSAT 在 求解 精度 方面 排名 第 一 ,并 且 与 CCEHC、 4 ”结束 语 
HS-Greedy 结果 相差 不 大 ， 但 远 高 于 了 SA、WalkSAT、 i 
CCLS2akms、Optirise6-in 几 种 算法 。 说 明 WWP+WalkSAT 算 本 文 设计 基于 变 元 权 值 计算 的 警示 传播 算法 WWP1+ 


+4 口 了 
日 口 
由 | 
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WalkSAT， 对 其 算法 进行 了 精度 分 析 ， 实 验 表明 该 算法 对 车辆 调度 、 机 器 人 路 径 规划 、 资 源 配 置 和 网 络 社交 分 析 等 
易 解 Max-3-SAT 问题 还 是 难 解 Max-3-SAT 问题 都 具有 较 方面 有 广阔 的 应 用 前 景 。 但 是 由 于 使 用 了 基于 变 元 权 值 计 
好 的 优势 ,而 且 也 打破 了 信息 传播 算法 求解 的 局 限 性 ， 其 算 的 警示 传播 算法 增加 了 整个 计算 的 复杂 度 因而 消耗 了 
结果 对 于 以 后 的 理论 研究 提供 了 很 大 的 参考 价值 。 与 此 同 时 间 , 后 续 的 研究 将 会 在 算法 的 变 元 剪 枝 策略 以 及 其 他 减 
时 该 算法 也 可 用 于 解决 不 同 领域 的 实际 问题 , 在 组 合拍 卖 、 少时 间 效 率 的 方法 入 手 研究 。 


表 2 SA、WalkSAT、WWP+WalkSAT、CCEHC、CCLS2akms、Optirise6-in 和 HS-Greedy 在 易 解 实例 中 实验 对 比 
Tab.2 Comparison of SA、walksat、 WWP+walksat、 CCEHC、CCLS2akms、Optirise6-in and HS-Greedyin easily solved instances 


Precision(The number of satisfy clauses / clauses) 


Instance Number 
SA WalkSAT WWP+WalkSAT CCLS2akms CCEHC Optirise6-in HS-Greedy 

s3v70c700 10 95% 92% 96% 93% 97% 86% 96% 
s3v70c800 10 95% 93% 96% 93% 96% 87% 96% 
s3v70c900 10 95% 92% 96% 93% 95% 86% 95% 
s3v70c1000 10 95% 93% 95% 93% 96% 88% 94% 
s3v70c1100 10 94% 93% 95% 93% 94% 88% 94% 
s3v90c700 10 96% 94% 97% 59% 97% 86% 96% 
s3v90c800 10 96% 94% 96% 59% 97% 86% 96% 
s3v90c900 10 95% 93% 96% 59% 96% 87% 95% 
s3v90c1000 10 95% 93% 96% 59% 96% 88% 96% 
s3v90c1100 10 95% 93% 96% 59% 96% 90% 96% 
s3v110c700 10 97% 95% 97% 60% 97% 86% 97% 
s3v110c800 10 97% 95% 99% 60% 97% 86% 97% 
s3v110c900 10 97% 95% 98% 60% 97% 88% 97% 
s3v1l10c1000 10 97% 95% 97% 60% 96% 88% 96% 
s3v110c1000 10 97% 95% 98% 60% 97% 87% 97% 


表 3 SA、WalkSAT、WWP+WalkSAT、CCEHC、CCLS2akms、Optirise6-in 和 HS-Greedy 在 难 解 实 例 中 实验 对 比 
Tab.3 Comparison of SA、walksat、 WWP+walksat、CCEHC、CCLS2akms、Optirise6-in and HS-Greedyin difficultly solved instances 


The number of unsatisfy clauses 


ee SA WalkSAT WWP+WalkSAT CCEHC CCLS2akms Optirise6-in HS-Greedy 
HG-3SAT-V250-C1000-1 23 64 9 5 N/A 138 12 
HG-3SAT-V250-C1000-2 25 66 10 6 N/A 122 13 
HG-3SAT-V250-C1000-3 27 64 9 6 N/A 129 10 
HG-3SAT-V250-C1000-4 28 61 8 5 N/A 115 12 
HG-3SAT-V250-C1000-5 24 69 10 6 N/A 110 11 
HG-3SAT-V250-C1000-6 27 60 4 6 N/A 128 14 
HG-3SAT-V250-C1000-7 23 66 6 6 N/A 120 14 
HG-3SAT-V250-C1000-8 26 64 10 3 N/A 123 13 
HG-3SAT-V250-C1000-9 22 65 9 7 N/A 115 12 
HG-3SAT-V250-C1000-10 27 68 5 6 N/A 114 13 
HG-3SAT-V300-C1200-1 29 68 13 6 N/A 156 18 
HG-3SAT-V300-C1200-2 28 67 10 6 N/A 155 20 
HG-3SAT-V300-C1200-3 25 70 13 y N/A 151 15 
HG-3SAT-V300-C1200-4 31 72 4 6 N/A 134 16 
HG-3SAT-V300-C1200-5 25 65 4 6 N/A 126 15 
HG-3SAT-V300-C1200-6 27 71 3 这 N/A 151 14 
HG-3SAT-V300-C1200-7 25 68 6 6 N/A 141 16 
HG-3SAT-V300-C1200-8 28 60 5 gE N/A 133 15 
HG-3SAT-V300-C1200-9 31 72 5 7 N/A 149 17 
HG-3SAT-V300-C1200-10 26 69 汉 6 N/A 155 18 
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